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Abstract 

An analytically solvable model is proposed exhibiting 1/f spectrum in any de- 
sirably wide range of frequency (but excluding the point / = 0). The model 
consists of pulses whose recurrence times obey an autoregressive process with 
very small damping. 
PACS: 05.40.+j, 02.50.-r, 72.70.+m 
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1. Introduction 

The puzzle of the origin and omnipresence of 1/f noise - also known as 'flicker' or 'pink' 
noise - is one of the oldest unsolved problem of the contemporary physics. Since the first 
observation of flicker noise in the currents of electron tubes more than 70 years ago by Johnson 
|I[ , fluctuations of signals and physical variables exhibiting behavior characterized by a power 
spectral density S (f) diverging at low frequencies like l/f s (5 ~ 1) have been discovered in 
large diversity of uncorrelated systems. We can mention here processes in condensed matter, 
traffic flow, quasar emissions and music, biological, evolution and artificial systems and even 
human cognition (see and references herein). 

1/f noise is an intermediate between the well understood white noise with no correlation in 
time and the random walk (Brownian motion) noise with no correlation between increments. 
The widespread occurrence of signals exhibiting power spectral density with 1/f behavior sug- 
gests that a general mathematical explanation of such an effect might exist. However, except 
for some formal mathematical descriptions like "fractional Brownian motion" or half- integral 
of a white noise signal || no generally recognized physical explanation of the ubiquity of 1/f 
noise is still proposed. Models of 1/f noise in some physical systems are usually specialized (see 
and references herein) and they do not explain the omnipresence of the processes with 
1/f 5 spectrum 

Note also some mathematical algorithms and models of the generation of the processes with 
1/f noise []10|-|i2|. These models can not, as a rule, be solved analytically and they do not 



reveal the origin as well as the necessary and sufficient conditions for the appearance of 1/f 
type fluctuations. 

History of the progress in different areas of physics points to the crucial influence of simple 
models on the understanding of new phenomena. We note here only the decisive influence of 
the Lorenz model as well as the logistic and standard (Chirikov) maps for understanding of 
the deterministic chaos and the quantum kicked rotor for revealing the quantum localization 
of classical chaos. 

It is the purpose of this letter to present and analyze the simplest analytically solvable 
model of 1/f noise which can be relevant for the understanding of the origin, main properties 
and parameter dependencies of flicker noise. Our model is a result of the search for necessary 
and sufficient conditions for the appearance of 1/f fluctuations in simple systems affected by 
random external perturbations, originated from the observation of a transition from chaotic to 



nonchaotic behavior in an ensemble of randomly driven systems JL3| , initiated in Ref . |L4| and 
further developed in Ref. |T5| . 

In the model there are analyzed currents or signals represented as sequences of random but 
correlated pulses whose recurrence times (intervals between transit times of pulses) obey an 
autoregressive process with small damping. It is shown that for small average recurrence time 
and very small damping, random increments of the recurrence times lead to 1/f behavior of 
the power spectrum of the signal or current in wide range of frequency, however, analytical at 
/ 

2. Model and solution 

Let us consider a point process when the intensity of some signal consisting from a sequence 
of pulses (elementary events) or current of particles through some Poincare section L m may be 
expressed as 

i{t) = Y, a5 (t-t k ). (i) 

k 

Here 5 (t) is the Dirac delta function, {tk} is a sequence of transit times tk at which the particles 
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or pulses cross the section L m and a is a contribution to the signal or current of one pulse or 
particle when it crosses the section L rn . 

The power spectral density of the current (1) is 



S {f)= hm l 2 4r 

^> T-,oc \ T 



j2 e -&*ftk 

h—h ■ 




J- \ (2) 



where A (k; q) = tfc+q— is the difference of transit times tk+ q and t^, T is the whole observation 
time interval, /c min and /c max are minimal and maximal values of index k in the interval of 
observation and the brackets (...) denote the averaging over realizations of the process. 

Let us analyze a process whose recurrence times Tk = tk — tk-i follow an autoregressive 
AR(1) process with offset f > 0, regression coefficient a = 1 — 7 and noise variance a 2 . So, 
defining by Ok — Tk — f a deviation of the recurrence time Tk from the steady state value f, we 
have autoregressive equations for the deviations 6k 

6 k = a&k-i + cre k (3a) 

and for the recurrence times Tk 

T k = T k -i - 7 {Tk-i -f) + otk- (36) 

Here {ek} denotes a sequence of uncorrected normally distributed random variables with zero 
expectation and unit variance (the white noise source) and a is the standard deviation of the 
white noise. Note that the coefficient 7 has a sense of damping (the relaxation rate of the 
recurrence times to the average value f ) and we will consider only processes with 7 <^ 1 or 
7 = 0. 

The recurrence equation for the transit times tk is 

tk = tk-i+ T k (4) 

where the recurrence times Tk is defined by Eq. (3b). 

The simplest interpretation of our model corresponds to one particle moving along some 
orbit. The period of this motion fluctuates (due to external random perturbations of the 
system's parameters) about some average value f. Some generalizations and extensions of the 
model and its interpretation will be discussed below. 

From Eqs. (3a) and (3b) follows explicit expressions for the deviation 9k and the recurrence 
time Tk, 

k 

9k = 9 a k + crY^ a k ~ j Ej, (5a) 

3=1 

k 

T k = f + (T -f)a k + aJ2 a fe ~% (56) 

j=i 

with 9 and r being the initial values of the deviation and of the recurrence time, respectively. 
The variance of the recurrence time Tk is 

al{k) = (Tl)-(Tk) 2 = ^(l-** k )/(l-**) (6) 
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After some algebra we can easily obtain from Eqs. (4) and (5b) an explicit expression for 
the transit times t k , 



t k = t + J2rj=t + kf + - (to - f) (l - a k ) + -J2 (l " e h 

7 7 i=i 



(7) 



with to being the initial time. 

The power spectral density of the current according to Eq. (2) depends on the statistics of 
the transit times difference A (k; q) which according to Eq. (7) is 



A (k; q) = t k+q -t k = qf + -(r - f) (1 - a q ) a k+1 

7 



+ - 



7 



k+q 



q>0. 



(8) 



(l-a«)J2a k+1 - l e l+ ]T (l - a k+q+1 ~ l ) e t 

1=1 l=k+l 

Random variable A (k; q) is a sum of two regular terms and k + q uncorrelated Gaussian 
random variables with zero expectations and variances (^j (1 — a q ) 2 a 2( - k+1 ~ 1 ^ for I — 1, 2, . . . , k 

and (j^J (l — a k+q+1 ~ 1 ^) for I = k + 1, k + 2, . . . , k + q, respectively. Therefore, A (k; q) is a 
normally distributed random variable with the expectation 



// A (k; q) ee (A (k; q)) = qf + - (r - r) (1 - a") a k+l 

7 



(9) 



and the variance a\ (k; q) = {A (k; q) j — (A (k; q)) , which equals the sum of the variances of 
the components, 

3=1 3=1 



a\ (fc; q) 



a 



J, 



Q- 



2a (1-ofl) a 2k+2 (l-a q Y 



1 -« 2 



1 -« 5 



q > 0. 



(10) 



Here new summation indexes j — k + 1 — I and j = /c + g + l — / of first and second sums in 
(8), respectively, have been introduced. 

Note that from the definition of the difference of transit times it follows the symmetry 
relations 



A (k; -q) — -A (k -q;q), H (k; -q) = -fi (k - q; q) , o\ (k; -q) = a\ (k - q; q) . (11) 

At k ^> 7" 1 or after averaging over r from the distribution with the expectation f and 
variance a 2 , ee a 2 (oo) = a 2 / (1 — a 2 ) ~ o" 2 /27 according to Eq. (6), expressions (7) and (8) 
generate a stationary time series. The expectation and the variance of the difference A (k; q) 
of transit times of the stationary time series are 



l_i A (q) ee (i_i A (oo;q)} = qr, 

2 



a 



Q- 



2a (1 - a q ) 
I -a 2 



(12) 
(13) 
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The power spectral density of the current according to Eq. (2) is 

2a 2 



S ( f) = lim ^- Y { e 

k,q 



i2TrfA(k;q) 



2a 2 

k,q 



(14) 



where XA(k-,q) (2vr/) is the characteristic function of the distribution of the transit times differ- 
ence A (k; q). For the normal distribution of A (k; q) the characteristic function takes the form 
XA(k-,q) (27r/) = exp [i2nf[iA (k; q) — 2n 2 f 2 <j\ (k; q)} and the power spectral density equals 



2a 



S (/) = lim — exp z27r//i A (k; q) - 2n 2 f 2 a 2 A (k; q) 

1 -^oo ± ' * 



(15) 



k,q 



For q <C 7 1 we have from Eqs. (9) and (10) expansions of the expectation /za (k; q) and of 
the variance a\ (k; q) in powers of 7 g <C 1 



(7 A (k; q) = 



fi A (k; q) =r (k) q, r (k) = r + (r - r) a 
'2(l- a 2k 



fc+i 



<7 



1 - a 2 



+ a 



2 k 



q + [a 



2k 



9* +3? 



The leading term of the expansion (17) may be written as 

a\ (k; q) = a 2 (k) q 2 , 7 g < 1 



(16) 
(17) 

(18) 



where the variance a 2 (k) of the recurrence time Tk is defined by Eq. (6). Therefore, Eq. (15) 
takes the form 

2a 2 



S (f) = lim — J2 exp i2nfr (k) q - 2n 2 f 2 a 2 T (k) q' 
1 -^00 ± * 



(19) 



k,q 



Here the mean recurrence time r (k) is defined in Eq. (16). 



Eq. (19) is valid if 27r 2 /V 2 (k) q 2 



l> 1, i.e., for f > fi= 7/7T(J r (k). When / < f T 



[2nr (k)] 1 and / < f 2 — [ncr T (k)] 1 we can replace the summation over q in Eq. (19) by the 
integration. The integration yields to the 1/f spectrum 



/ V 7T T^oo T ^ a T (k) 



exp 



2a 2 (k) 



, f\ < f < min{/ 2 ,/ T } 



(20) 



Eq. (20) is valid for the stationary as well as for the non- stationary process with slowly 
changing of the expectation and variance of the recurrence time. 

At k ^> 7 _1 or after averaging over To from the distribution with the expectation f and 
variance a 2 = a 2 / (1 — a 2 ) ~ o- 2 /2 7 we have the stationary process: the expectation and 
the variance of the recurrence time do not depend on the parameter k, i.e., r (k) = f and 
a T (k) = a T . For the stationary process Eq. (20) takes the form 



s(f) = i 2 J- 



(21) 



Here / = lim a (/c max — &mm + 1) /T — a/r is the average current and an is a dimensionless 

T—>oo 

constant (the Hooge parameter) 



Therefore, the power of 1/f noise except of the squared average current strongly depends 
on the ratio of the average recurrence time to the standard deviation of the recurrence time. 
3. Discussion and generalizations 

The point process containing only one relaxation time 7 _1 can for sufficiently small damp- 
ing 7 and average recurrence time f <C G I \fl (with o being the standard deviation of the 
white noise source) produce an exact 1/f-like spectrum in wide range of frequency (/i,/2), 
with /2//1 — 7 -1 - Furthermore, due to the contribution to the transit times of the large 
number of the random variables E\ (I = 1,2, ...k), our model represents a 'long-memory' ran- 
dom process. As a result of the nonzero relaxation rate (7 7^ 0) and, consequently, due to 
the finite variance, a 2 = a 2 /2 / y, of the recurrence time the model is free from the unphysical 
divergency of the spectrum at / — > 0. So, using an expansion of expression (13) at 3> 1 
, cta(<?) = h) 2 <!, we obtain from Eq. (15) the spectrum density S(f) = I 2 (2a 2 /f / y 2 ) for 



/ <C min{/i,/o — f^/ncr 2 }. This is in agreement with the statement ]Tj| that the power 
spectrum of any pulse sequence is white at low enough frequencies. 

This simple, consistent and exactly solvable model can easily be generalized in different 
directions: for large number of particles moving in similar orbits with coherent (identical for all 
particles) or independent (uncorrelated for different particles) fluctuations of the periods, for 
non-Gaussian or continuous perturbations of the systems' parameters, for nonlinear relaxation 
and for spatially extended systems. So, when an ensemble of N particles moves on closed 
orbits and the period of each particle fluctuates independently (due to the perturbations by 
uncorrelated sequences of random variables {£%}, different for each particle v ) the power spectral 
density of the collective current / of all particles can be calculated by the above method too 
and is expressed as the Hooge formula 0] 



S(f) 



f2 a H 



Nf 



(23) 



The model may be used for evaluation of the power spectral density of the non-stationary 
process as well. So, at k <C 7 -1 or for 7 = we have a process with the constant averaged 
recurrence time r (k) = r , from Eq. (16), and linearly increasing variance of the recurrence 
time a 2 (k) = a 2 k, according to Eq. (6), i.e., a process similar to the Brownian motion without 
relaxation. For k 3> \q\ it is valid expansion (18) and the power spectral density of such process 
for finite observation time interval T may be evaluated according to Eq. (20) 



S(f) 




E 



Vk 



cxp 



2a 2 k 



h < / < / 2 



(24) 



Here 



h 



ak n 



ku 



-1 



h 



1x0 



-1 



(25) 



The process with random increments of the recurrence time and without relaxation is, 
however, very unstable and strongly depending on the realization. Really, the averaged number, 
(&max — &min) , of transition times tk for the given observation time interval T is (fc max — ^min) = 
T/tq but the standard deviation <jt of the time interval for given fc min and fc max according to Eq. 

> 3 (ro/of 



(10) equals a T = a (k max - k min ) yj (k max + 2k min ) /3. Therefore, for (fc max + 2fc min ) 
the standard deviation exceeds the expectation value (T) = (/c max — k m[n ) r of the time interval. 
The power spectrum of any realization is, however, of 1/f type for large frequency interval 
C/1,/2), with f 2 /ft ~ Jk 

max^min- 
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It should be noticed in conclusion that in many cases the intensity of signals or currents 
can be expressed in the form (1). This expression represent exactly the flow of identical point 
objects. More generally, in Eq. (1) instead of the Dirac delta function one should introduce 
time dependent pulse amplitudes A k (t — t k ). However, the low frequency power spectral density 
depends weakly on the shapes of the pulses |IB| , while fluctuations of the pulses amplitudes 
result, as a rule, in white or Lorentzian, but not 1/f, noise. The model (2)-(5) in such cases 
represents fluctuations of the averaged period r k between the subsequent transition times of the 
pulses. Therefore, the model may be easily generalized and applied for the explanation of 1/f 
noise in different systems. Furthermore, it reveals the possible origin of 1/f noise, i.e. random 
increments of the time intervals between the pulses or elementary events. 

Summarizing, a simple analytically solvable model of 1/f noise is presented and analyzed. 
The model reveals main features and parameter dependences of the power spectral density of 
the noise. The model and its generalizations may essentially influence the understanding of the 
origin and main properties of the flicker noise. 
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